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Cn 1 1 INTRODUCTION 



ABSTRACT 

There are only a few ways to constrain the Era of Reionization and the properties 
of high redshift (z > 6) stars through observations. Here, we discuss one of these 
observables - the spectrum of the Near Infrared Background - and how it is po- 
tentially affected by the transition from Population III to Population II stars. The 
stronger Lyman-a emission expected from massive Population III stars could result in 
a 'bump' in the spectrum of the Near Infrared Background (referred to in this work 
as the Lyman-a bump). The strength and shape of this bump can reveal properties 
of Population III stars. The Lyman-a bump is predicted to be higher if Population 
III stars are more massive and present at lower redshifts. The shape of the bump 
is governed by the star formation rate and the time it takes Population III stars to 
transition to Population II stars. If Population III stars are indeed massive, a bump is 
predicted as long as Population III stars exist at z < 15, even if their star formation 
rate is as low as 10~ 7 M Q yr~ 1 Mpc -3 . This means that there may be some observa- 
tional signature in the Near Infrared Background of small pockets of metal-free gas 
forming Population III stars at z ~ 6, even if they are quite rare. 
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The first generation of stars, known as Population III, were 
truly metal free. These stars have never been observed, so 
their properties must be surmised using indirect or the- 
oretical means. It is possible that these true first gen- 
eration of stars were on the order of hundreds of sola r 
masses l|Abel et al.ll2002l : iBromm et alJliool lGloverll20o3) . 
although some s uggest these star s were less massive with 
a broader IMF (IClark et all 120081 l201ll ; IStacv et all l20ld : 
iGreif et al.ll201 J : Hosokawa et alj|201ll h As star formation 
progressed and these stars died, the intergalactic medium 
was enriched with metals. After several generations of en- 
richment, the amount of metals present would be signif- 
icant enough to act as an efficient coolant, causing the 
properties of subsequent generations of stars to be funda- 
mentally different. These Population II stars, forming in 
clouds that can undergo metal and dust cooling, were able 
to fragment into smaller masses, likely causing the masses of 
these stars to be less than their Population III progenitors 
(IBromm et al.ll2002l ; ISchneider et al.ll2002l : IBromm fc Loebl 
l2003h . Currently, it is unknown when this transition from 
Population III to Population II stars occurred and how ex- 
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tended in time it was. However, even though the majority 
of star formation transitioned to Population II, it is possi- 
ble that Population III formation persisted even to z ~ 6 in 
pock ets of pristine gas |Tornatore et alj 120071 ; iTrenti et al.l 
120091) . 

Observing this era of high-redshift star formation is 
very challenging. High redshift galaxy surveys are now rou- 
tinely discovering galaxies at very high redshifts (z > 6). 
However, these detections are currently limited to only the 
brightest objects common enough to appear in the survey 
field, missing out on the bulk of numerous, smaller galax- 
ies. These galaxies below the detection limit could be re- 
sponsible for a significant amount of star formation, and 
could be the primary drivers of reionization jBouwens et al.l 
l20ld : iRobertson et al.ll20ld ; iFernandez fc Shulfeoul h Even 
though these galaxies are below the limiting magnitude 
of current surveys, it may be possible to observe their 
redshifted cumulative spectra, which would be present in 
any background emission in the infrared. Therefore, under- 
standing what component of the Near Infrared Background 
(NIRB) is due to a high redshift component could lead to 
important constraints on these early stars and galaxies. 

Actually measuring the contribution of these high 
redshift stars to the NIRB is no easy task, due to great 
uncertainties in foreground subtraction. In order to accu- 
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rately establish which portion of the spectrum is actually 
due to stars during the era of reionization, low redshift 
galaxies must be correctly subtracted. Zodiacal light, 
which is very difficult to model, must also be taken into 
account. Yet, despite the difficulties, many have undertaken 
the challenging observation to measure this excess of the 
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2005T). These observations remain controversial, with some 



claiming t hat the NIRB is r e solved - either in lower redshift 
galaxies dThompson et all l2007al lbh or intrahalo stars 
l|Coorav et al.ll2012bh . ~ 



However 
reionize the 
Reionization would 
NIRB, even if it is 



because energetic photons are needed to 
universe, any signal from the Epoch of 
be redshifted and present in the 
not detectable with current instru- 



ments. In this spirit, many have sought to model the 
contribution of stars and galaxies responsib le for reion- 
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vations, such as from CIBER, AKARI, and CANDLES, 
better constraints on any observable contribution of stars 
from the Epoch of Reionization to the NIRB should be 
better established. 

In this work, we explore how the transition from Popu- 
lation III to Population II stars is imprinted on the spectrum 
of the NIRB. Our stellar models, in addition to describing 
the star formation rate and how we deal with the transition 
from Population III to Population II stars, are described in 
Sj2] Our formalism for obtaining the spectrum of the NIRB 
is detailed in §3] The effect of various parameters on the 
NIRB spectra, such as the mass of stars, the dependence 
the star formation rate, the contribution of the IGM, and 
parameters of the Population III to Population II transition, 
are detailed in <J4] We conclude in Sj5] 



2 STELLAR POPULATION MODELS 

In order to understand the spectra of high redshift stars, 
we need to establish stellar properties, such as mass and 
metallicity. This metallicity will increase as time goes on, 
so a model for the transition from metal-free to metal-poor 
stars must be established. We also must model how many 
stars are forming and how this changes with redshift. 



2.1 Mass and Metallicity of High Redshift Stars 

In order to model the various stellar populations, we estab- 
lish some limiting cases for both the mass and metallicity 
of our stars. For the metallicity of the stars, we use either 



Population III stars with no metals, or Population II metal- 
oor s tars (with Z = 1/50 Z ), as in IFernandez fc Komatsul 
2006). Since the mass of Population III stars is still largely 

unknown, we use three possibilities for the initial mass spec- 

tra. The first two are represented by a Larson mass spectra 

(|Larsonlll998T l: 



f(m) oc m 1 ( 1 + — 
' m, 



(1) 



For very massive stars, we set the mass limits to mi = 
0.1 M Q , m 2 = 500 M , and m c = 250 M Q , while for less 
massive stars, we use mi = 0.1 Mq, 7712 = 150 Mq, and 
m c = 10 Mq. For the case with the least massive Popu- 
lation III stars, along with Populatio n II stars, we use a 
Salpeter mass function i|Salpeterlll955r i: 



,1 \ -2.35 
/(m) oc m , 

with mass limits of mi 



(2) 



0.1 M Q and m 2 = 150 M e 



2.2 The Star Formation Rate 

The star formation rate at high redshifts is very hard 
to probe observationally . There are some result s deduced 
from gamma-ray bursts dRobertson fc Ellis!l2012r) and from 
high redshift galaxy surveys llHopkins fc Beacoml 120061 : 
iMannucci et al.ll2007l : iBouwens et al.ll2008l ) up to a redshift 
of z ~ 9. However, these results may not probe the popula- 
tion as a whole (including the smallest of galaxies), and the 
star formation rate at even higher redshifts is still uncon- 
strained. 

Knowledge of the star formation rate, then, is mostly 
obtained through theoretical means. In Fig. [1] we show 
the star formation rate based off analytical models from 
iTrenti fc Stiavellil ([2009) , who modeled star formation based 
on the collapse fraction of haloes, along with gas cooling, ra- 
diative feedback, and metal enrichment. The data was taken 
from their standard model, which contains only one Popu- 
lation III star per halo (purple dashed-triple dotted lines), 
or multiple Population III stars per halo (green long dashed 
lines). They calculated the values of the star formation rate 
of Population II stars (top lines for each colour) and Popu- 
lation III stars in either mini-haloes or more massive haloes 
with T V i r ^ 10 4 AT (bottom lines of each colour). We also 
show result s from some recent simulations. The cyan dot 
dashed line (|johnson et al.ll2013h shows the result of cosmo- 
logical simulations that contain metal e nrichment and stellar 
radiation fields. The simulations from () Alvarez et al.l[2012t ) 
(the red short dashed line) show the star formation rate for 
Population II stars only from simulations which take into ac- 
count suppression on small galaxies within ionized regions. 

In order to model the star formation rate p(z), we ap- 
proximate a fit as: 



p(z) = 10 



!/0+!/l z + I/2Z 



(3) 



where 1/0, yi, and 1/2 are constants. This model is not meant 
to be physical, but instead it approximates the star forma- 
tion rate to be consistent with results from analytical formu- 
lations and simulations. We set yi — —0.03, and change 2/2 to 
be either -0.002, -0.004, or -0.006. This effectively changes 
the steepness of the slope of the star formation rate, where 
more negative values of j/2 lead to less star formation at high 
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Figure 1. Our star format ion rate models, compa red with re- 
sults from analytical models llTrenti fc Stia vclli 2009) and simula- 
tions l lAlvarez et al J2012l ; |johnson et al.l2013h . The parametrized 
models are the models we use, based on equation \3\ which are 
shown as the solid lines for y 2 = -0.002, -0.004, or -0.006. The 
pink dotted line shows our second model for the star formation 
rate when a constant value of /* = 0.003 and a star formation 
time-scale of 11.5 Myr is assumed. 



io-' r 




z 



Figure 2. The total star formation rate of Population III and 
Population II stars with j/2 = —0.004 (solid black line) and Pop- 
ulation III stars only (coloured lines) for various values of the 
redshift of transition from Population III to Population II stars 
(zt) and the transition width (<r p ). 



redshift. The value of yo is then set so that the star formation 
rate at z — 6 is 0.02 MQyr _1 Mpc~ 3 , approximately con- 
sistent with current observational constraints. These models 
are shown as the solid black lines, also in Fig. [T] As a final 
model for the star formation rate, we also assume a constant 
star formation efficiency over redshift, which has often been 
assumed in the literature. This is shown as the dotted pink 
line. (This model will be discussed more in Section T4.1 1) . 

2.3 The Transition from Population III to 
Population II Stars 

We assume that the transition of Population III to Popula- 
tion II stars is described by the formula: 

U = \[l + erf{^)] (4) 

|Coorav et al.l l2012al ). This formula quantifies the (likely 
very in-homogeneous) transition with only a few free param- 
eters. We then define the star formation rate of Population 
III stars as f p p(z) and the star formation rate of Population 
II stars as (1 — fp)p(z). These are parametrized by the red- 
shift of transition (zt), and the width of the transition (a p ). 
Both zt and a p can be adjusted to simulate various enrich- 
ment histories. We show our total star formation rate (given 
by the solid black line in Fig. [2} for yi — —0.0004, along with 
the star formation rate of Population III stars only (given 
as the coloured lines) representing various values of zt and 
a p . 

2.4 Constraints from Reionization 

Properties of the stars themselves, such as the number of 
ionizing photons produced per second per star, Q(H), and 
the lifeti me of th e star, r*, are taken from stellar mod- 
els from ISchaererl ((2002) , or from the fitting formulas in 



iFernandez fc Komatsul (|200cf ) based on the stellar models 
in lLeieune fc Schaererl (|200ll ). When combining these mod- 
els with the mass spectrum f(m) defined in section [2~T1 we 
can compute the rate of ionizing photons produced per co- 
moving volume. The formula used depends on the mass of 
the star. Stars that are massive (above a mass mt(z), whose 
lifetimes are shorter than the star formation time-scale, here 
assumed to be the age of the universe) have a photon pro- 
duction rate based off the star formation rate p(z): 



r _ C: {z) Q{H)r{m)f{m)dm _ 

$21 dm f( m ) m 



(5) 



Less massive stars (denoted by the subscript s), whose life- 
times are longer than the age of the universe, do not emit all 
of their ionizing photons over the time interval we are con- 
cerned with, thus, the number of ionizing photons per stellar 
mass is computed instead by using the comoving mass den- 
sity of stars (p(z)): 



T s = 



Q Q(H)f(m)dm 
XT,' dmf{m)m 



P*(z) 



(6) 



We assume that the ratio of ionizing photons that es- 
cape into the IGM by the end of reionization to baryons 
T) is 2, which is know n as photon-starved reionization 
Bolton fc Haehnelti r2007). We can then solve for the escape 
fraction over the interval zi < z < 22: 



T 



n H (l + z)H(z) 
Tiii{z) + Tii(z) 



dz. 



Therefore, the number density of hydrogen (n«) is balanced 
by the number of photons produced by each population of 
stars, where Tin or Vn is the sum of Ti and F a for either 
Population III or Population II stars. 
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3 THE INTENSITY OF THE NEAR 
INFRARED BACKGROUND 

Now that our stellar models are established, we can com- 
pute the portion of the NIRB expected due to high redshift 
stars within our redshift range (6 < z < 30). This emis- 
sion would result from a combination of various stellar and 
nebular processes. The intensity of the NIRB is given as: 

c [ dzp([l + z]v, z) 

47T 



(8) 



H(z)(l + z) ' 

jPeacockl Il999h . In order to compute the component 
of the intensity of the Near Infrared Background from 
high redshift stars, we follow the formalism presented in 
iFernandez fc Komatsul {2006). Again, in order to compute 
the emissivity p([l + z\u, z), we must divide our calculation 
into parts, depending on the lifetime of the star. If the life- 
time of the star is shorter than the star formation time-scale, 
the emissivity is given by the luminosity L a (m) of each stel- 
lar or nebular component a, integrated over a mass spec- 
trum of stars f(m), and is weighted by the lifetime of the 
star r(m). This is normalized by to*, or the mean stellar 
mass of the mass spectrum: 



p(v, z) = p*(z)c / dm mf(m) 



L"(m)r(m) 



(9) 



and if the lifetime of the star is longer than the star forma- 
tion time-scale, the emissivity is given by: 



p(v,z) = p,(z)c 2 — 



dm mf(m) 



L a v {m) 



(10) 



(See appendix A of IFernandez fc Komatsul (|2006l )V Emis- 
sion will originate from the star itself (in the form of a 
stellar black-body), as well as nebular emission, namely, 
free-free, free-bound, and two-photon emission, along with 
emission lines, the primary one being the Lyman-alpha line. 
The luminosity of each comp onent L™(m) is modeled as in 
IFernandez fc Komatsul {2006). For our initial cases, we as- 
sume that only emission from the haloes contribute to the 
spectra of the NIRB, so the luminosity of the nebular com- 
ponents will be multiplied by the factor (1 — fesc)- 



4 THE SPECTRAL FEATURES OF THE NIRB 

For our reference model, we allow star formation to begin 
with only Population III stars with a Larson mass spectrum 
and m c of 250 (see equation [TJ. At zt = 10, star formation 
transitions to Population II stars via equation|H with a tran- 
sition width of <t p = 2. The star formation rate is assumed 
to follow equation [3] with 1/2 = —0.004, and / esc is set to 
be consistent with photon-starved reionization (in this case, 
fesc = 0.16). These assumptions (the masses of the stars, 
details of the transition from Population III to Population 
II stars, and the star formation rate) will be adjusted in the 
following sections. 



4.1 Moving from Population III to Population II 
Star Formation 

It is unknown when the majority of star formation transi- 
tioned from Population III to Population II stars. In Fig. 



[3] we allow this redshift of transition (zt) to be either 8, 
10, 12, or 15. For our reference case, Population III stars 
are very massive and will have very strong nebular features, 
such as the Lyman-a line. Because of this, it is possible 
that the overall spectral intensity for a more distant popu- 
lation of Population III stars could have a higher intensity 
than closer, Population II stars. Because of this, the spec- 
trum of the NIRB could contain a 'bump' from the Lyman-a 
line from Population III stars at higher redshift than closer, 
Population II stars. (A simila r prominent Lyman-a feature 
is seen bv iDwek et al. I (|2005l 'l. who assume all stars form- 
ing at high redshift are Population III stars.) If zt is larger, 
the Lyman-a' bump is less pronounced since these Popu- 
lation III stars are more distant, and the overall spectrum 
of the NIRB becomes increasingly dominated by Population 
II stars with a less prominent Lyman-a line. As z t rises, 
the spectra of the NIRB becomes increasingly featureless at 
longer wavelengths. If zt > 15, no Lyman-a bump is seen 
for our given star formation rate and initial mass spectrum. 

This has interesting observational consequences. If there 
are some Population III stars at lower redshift (and even 
if the star formation rate of these stars is quite low) they 
could contribute a signature to the NIRB spectrum, causing 
it to deviate from a featureless spectrum. For example, our 
reference case has pin(z = 6) ~ 10 -5 M0yr -1 Mpc~ 3 , (in 
other words, the Population III star formation rate is 0.2 
per cent of the total star formation rate at z — 6). Our case 
where z t = 12 has pm(z = 6) ~ 10" 7 M©yr _:L Mpc~ 3 , or a 
Population III star formation rate of 0.001 per cent of the 
total at 2 = 6. These star formation rates are similar to rates 
predicted fro m pockets of metal free gas that may ex ist at 
low redshift {Tornatore et alj|2007l : iTrenti et alj|2009h . 

The transition from Population III to Population II 
stars is also described by the width of the transition, 
parametrized by a v . In Fig. U a p is adjusted, so as to cause 
the transition to be abrupt (<r p = 0.5) or extended in time 
(a p = 10). When a v — 10, there are still significant Popu- 
lation III stars forming at 2 ~ 6. Because of this, there is 
significant Lyman-a emission at shorter wavelengths from 
these massive stars. If a v is smaller, this low-redshift contri- 
bution from massive Population III Lyman-a emission falls, 
causing a drastic change in shape of the Lyman-a bump. 
If the transition is very abrupt, there can actually be more 
emission from more distant Population III stars than nearby 
Population II stars, creating a situation where there is more 
emission at longer wavelengths than shorter wavelengths. 

Also shown in Fig.|4]is the effect of changing the slope of 
the star formation rate as a function of redshift, described by 
our free parameter j/2 (in equation [3]), where more negative 
values of j/2 lead to a more rapidly decreasing value of the 
star formation rate. Changing this free parameter has only 
a slight effect on the spectra of the NIRB. This is because 
changing 7/2 has the greatest effect at high redshift (since 
equation [3] is normalized at z ~ 6). The stars at higher red- 
shift, since they are more distant, have a smaller contribu- 
tion on the overall shape of the spectra of the NIRB. There- 
fore, changing a v can change the shape of the Lyman— a 
bump. 

However, instead of assuming a star formation rate as 
a function of redshift via equation [3] we can instead assume 
the fraction of baryons that form into stars (/*) is constant 
with redshift. (This method was used in previous papers, 
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Figure 3. The spectra of the NIRB, adjusting the redshift where the majority of star formation transitions from Population III to 
Population II (zt). If zt is relatively low, a 'bump' in the spectra is seen from the end of the Population III epoch. The solid lines are 
the total contribution from z = 6 — 30, dot-dashed are z = 6 — 8, dashed are z = 8 — 10, dotted are z = 10 — 12, and dashcd-triplc dot 
are z = 12 - 18. 




5.0 3.0 2.0 1.0 0.7 0.4 

Observed Wovelength [microns] 



Figure 4. The effect of the spectra of the NIRB from changing a p , which describes the width of the transition from Population III to 
Population II stars. Also shown is the dependence on the star formation rate, described by the parameter j/2i which affects the steepness 
of the star formation rate as a function of z, or by assuming instead a set fraction of baryons (/*) that form into stars over time. 
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suc h as Fernandez et alJ (|201dt ). iFernandez et al.l (|2012bh . 
and lCoorav et all l|2012aT l). The star formation rate equiva- 
lent to this assumption is shown by the pink dotted line in 
Fig. □ for /» = O.OOfl derived by combining /„ with the 
collapse fraction of haloes, and as suming a time formatio n 
time-scale (here, 11.5 Myr, like in IFernandez et all l|20ld) ). 
The star formation rate as a function of redshift for this 
case is much steeper. Because of the lower level of the star 
formation rate at higher redshift, the Lyman-a bump is not 
present, because lower redshift Population II stars will domi- 
nate over any Population III contribution. The spectrum for 
this case is shown as the blue dashed-dotted line in Fig. U 
It is important to note that the assumption of a constant 
/, , as often assumed in the literature, predicts a lower star 
formation rate at high redshift in relation to lower redshifts 
as compared to our other models. 



0.100 



0.010 



0.001 
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4.2 The Masses of Population III Stars 

To model the possibility that Population III stars might not 
be as massive, we computed the spectrum of the NIRB when 
the value of m c for Population III stars is 10 Mq and the 
mass limits of the mass spectrum are mi = 0.1 Mq and 
— 150 Mq (see equation [TJ. In addition, we also calcu- 
late the spectra of the NIRB if Population III stars have a 
Salpeter slope with mi = 0.1 Mq and m,2 = 150 Mq. The 
results of this test, in comparison to our reference case with 
heavier Population III stars, are shown in Fig. [5] 

When Population III stars are more massive, there is 
a more pronounced bump from the Lyman-a line from the 
end of the Population III era. But, if Population III stars 
are no more massive than Population II stars, the spectra 
is smooth, and the Lyman-a bump is not seen at all. More 
massive Population III stars would emit a stronger Lyman- 
a emission line, therefore the presence of a Lyman-a bump 
would be the hallmark of a massive Population III era. 

4.3 The Contribution of the IGM and the Escape 
Fraction 

Nebular emission will originate from within the halo, or, 
if some photons escape the halo, from the IGM. However, 
because of the low density of the IGM, the IGM component 
of the NIRB might be quite low. Therefore, it is possible that 
the contribution of the IGM to the overall spectrum of the 
NIRB can be negligible (jNakamoto et al ]|200ll : ICoora7 et al.l 
l2012al ). 

Thus far, we have assumed that the IGM does not con- 
tribute at all to the overall intensity of the NIRB. In this 
case, nebular emission will only originate from the haloes 
themselves, so their luminosity must be multiplied by the 
factor (1 — fesc)- If, on the other hand, the recombination 
time within the IGM is very short, the IGM will contribute 
to the spectra of the NIRB. As a comparison, we show the 
spectra of the IGM if the entire IGM is contributing to the 
NIRB in Fig. [6] In this case, all photons contribute to the 
spectra of the NIRB, whether they are from the haloes or 



This value was chosen so that the star formation rate is ap- 
proximately consistent with our other models at z ~ 6. 



Figure 5. The spectrum of the NIRB when Population III 
stars have various stellar masses. In each case, the Population 

II stars have the same mass spectrum. These are combined with 
Population III stars with a heavy Larson mass spectrum (with 
m c = 250 Mq - dashed lines), Population III stars with a light 
Larson mass spectrum (with m c = 10 Mq - solid lines), and Pop- 
ulation III stars with a Salpeter mass spectrum (dotted lines). 

the IGM, so the nebular component of the emission is no 
longer multiplied by the factor (1 — / esc ). 

The spectra of the NIRB is very weakly dependent on 
whether or not the IGM is contributing to the overall emis- 
sion. This is because, for this case, the escape fraction of 
photons is only 0.16, set to be consistent with constraints 
from reionization. Therefore, the nebular component of the 
emission is only diminished by ~ 16 per cent when the IGM 
is not contributing at all to the emission in the NIRB. This 
is not enough to make a significant change in the spectrum 
of the NIRB. 

Yet, changing the value of the escape frac tion can affect 
the spectrum of the NIRB. lYue et ail (|2012T l did not see a 
Lyman-a bump in their predicted spectrum of the NIRB. 
This was because they allowed fesc to vary with redshift. 
According to their model, f eac (z = 5) = 0.05, and rises 
with redshift, reaching unity at z ~ 11. Because of this, any 
Lyman-a contribution from stars at higher redshift would 
be suppressed, and therefore, no Lyman-a bump from these 
stars would be seen. Therefore, the presence of a Lyman-a 
bump would also indicate that / e3C is low enough to allow a 
prominent Lyman-a feature at higher redshifts. 

It is also possible that the escape fraction for Population 

III stars is above that of Population II stars (especially if 
the Population III stars are significantly massive) . This case 
would further smooth out any effect of the Lyman-a bump, 
since the Lyman-a emission from higher redshift Population 
III stars in haloes would be diminished. 



5 CONCLUSIONS 

The spectra of the Near Infrared Background can give infor- 
mation on the end of the Population III era. This is mainly 
reflected in the presence and shape of the 'Lyman-a bump', 
whose shape reflects not only when Population III stars tran- 
sitioned to Population II stars, but also properties of these 
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Figure 6. The spectrum of the NIRB, with and without contri- 
bution from the IGM. 



stars as well. The presence of this bump is mostly dictated 
by the relative presence of massive Population III stars to 
less massive Population II stars. A larger bump corresponds 
to a situation where more massive Population III stars exist 
to lower redshifts. For our reference case, this bump in the 
spectrum is seen when zt < 15 (or, in other words, a case 
where the Population III star formation rate is greater than 
~ 1CP 10 times that of the total star formation rate). 

The shape of the Lyman-a bump also will reveal infor- 
mation on how quickly the Population III era ended. An ex- 
tended tail of Population III to low redshifts, which would 
result from a very long transition from Population III to 
Population II stars, would result in a strong Lyman-a peak. 
However, if the transition is abrupt, the Lyman-a bump will 
change shape, such that it is possible that there is more in- 
tensity at longer wavelengths than shorter wavelengths, a 
result of more emission from more massive, more distant 
Population III stars than nearer, less massive Population II 
stars. 

However, if no Lyman-a bump is seen at all, this indi- 
cates that the Population III era ended much earlier, with 
zt > 15, or that Population III stars are not very massive, 
being closer to mass to the mass spectrum we see today. On 
the other hand, it can also indicate that the escape fraction 
increases at high redshifts so that the Lyman-a emission 
from higher redshift Population III stars is suppressed, and 
that the emission from the IGM is negligible. 

Even though the mean NIRB is difficult to observe, the 
relative change in intensity resulting from a Lyman-a bump 
as a function of wavelength will show up as a deviation from 
a flat continuum. It may be possible to observe this devi- 
ation, independent of complete subtraction of foregrounds, 
so long as these foregrounds have a flat spectra. 

These results have interesting consequences. It is pos- 
sible that small pockets of metal-free gas persist until low 
redshifts, so that Population HI star f o rmation might sti ll 
persist at z ~ 6 (jTornatore et al.ll2007l ; iTrenti et alj|2009l ). 
Even small Population III star formation rates (p(z = 6) ~ 
4 x 1CT 5 or 2 x 10~ 7 M yr _1 Mpc~ 3 for our cases a t = 2 
and zt — 10 and 12 respectively), could actually lead to a 
signature on the spectrum of the Near Infrared Background. 



Of course, if these pockets of metal free gas are indeed rare 
on the sky, measurements must be carefully done to assure 
this effect is observed. 
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